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Abstract 

We study the dynamics of 3 point-vortices on the plane for a fluid 



> 

■ governed by Euler's equations, concentrating on the case when the 
I moment of inertia is zero. We prove that the only motions that lead to 

^SJ ' total collisions are self-similar and that there are no binary collisions. 

Also, we give a regularization of the reduced system around coUinear 

■ configurations (excluding binary collisions) which smoothes out the 

dynamics. 



(ii ! 1 Introduction 

> 

The dynamical equation for the vorticity field co on an ideal incompressible 
fluid with no viscosity is derived from Euler's equation and is given by 
■ Duj/Dt = uj ■ Vu. Here u is the velocity field, u; = V x u is the vorticity field 

^ . and D/Dt := d/dt+u-V is the material derivative. The context of this paper 

^ I is the situation when the vorticity field u; is highly localized around a finite 

number of points, in effect being described by a delta function. This is known 
as the A/^-point- vortex (or A'^-vortex, for brevity) problem. Each point-vortex 
is completely characterized by its position and its vortex-intensity. It is 
analogous to the iV-point-masses considered in celestial mechanics, except 
that, in contrast to the masses, the vortex-intensities can be positive or 
negative. We will concentrate on the case when the motion takes place on 
the plane. This corresponds to the physical situation of parallel vortex 
filaments on three-dimensional space. 
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The equations of motion for the A^-vortex problem on the plane, given 
by are derived in most standard texts on hydrodynamics, e.g. [Sj- It 
is a Hamiltonian system. Its solution for = 2 is a routine exercise. For 
A'^ > 4 it is a non-integrable system. The case = 3 is interesting because 
it is still integrable and it has been shown to be a fundamental component 
of the A^-vortex problem for large A^. 

The study of the 3 point-vortex problem was pioneered by W. Grobli in 
his 1877 dissertation^, who established the integrability of the system and 
obtained analytical solutions for certain special cases. Several decades later, 
J. L. Synge |3] studied the dynamics of the reduced system obtained by ig- 
noring rotations and dilations, giving a classification of the possible reduced 
dynamics based on the harmonic mean (the virial) of the vorticities. For 
this he used a geometric method based on using the so called trilinear coor- 
dinates and an argument based on the linearization around the equilateral 
configuration. Later, E. A. Novikov [£] gave an account of the solutions to 
the problem with three equal vorticities and, with Yu. B. Sedov studied 
self-similar collapse (defined below) and the dynamics near collapse when 
the harmonic mean of the vorticities is zero. H. Aref ^ gave a further 
classification of the possible motions according to the harmonic, geomet- 
ric and arithmetic mean of the vorticities, bringing together the methods 
in IHl and (S]. In another significant development, D. M. F. Chapman |^ 
studied the relation between conserved quantities and symmetries (transla- 
tional, rotational, temporal and scale transformations). Most of the results 
in these references are nicely summarized and discussed in §2.2], where 
some additional material can be found. 

This paper deals with the 3-vortex problem on the plane and has two 
main purposes. The first one is to give a negative answer to the question 
of whether there exist non-self-similar collisions. Surprisingly, this question 
has not been fully addressed in the literature. (Aref 's []P| description of the 
dynamics for zero moment of inertia does imply the impossibility of non- 
self-similar trajectories leading to collision, but his presentation is based 
on illustrative examples and not a definite proof.) As a by-product of our 
answer to this question we offer a simple description of the dynamics for 
the case when the moment of inertia is zero (which contains the self-similar 
collisions when the virial is null). This description necessarily overlaps with 
some results in the references mentioned above, but the method that we 
use, based on a projective blow-up, is different. The second objective of this 
paper is to describe the construction of a coordinate system that regularizes 

^Cf. 2 for a detailed account on Grobli's work. 
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the vector field of this system near collinear configurations. This is another 
sort of blow-up which duplicates the phase space of a vector field around 
singular points where a square root becomes zero. As part of this discussion 
we also show that for the three vortex problem the binary collisions are 
excluded from the dynamics. 

The system of N point-vortices on the plane can be described as follows. 
Let = + iya G C be the coordinates of the a-th point vortex with 
vorticity Ta, a = 1, . . . , N. Thus the configuration space is C" = M?^. The 
equations of motion are given by^ 



N 

This system is Hamiltonian, which can be seen as follows. Give C" the 
symplectic form O defined by 

n 
a=l 

where the E M, Fq, 7^ 0. Then Hamilton's equation ixn^ — takes 
the form 

OH 

F„i„ = -2i^ , a = l,...N . (2) 



Let ^ 

^ ~ ToF/j In lap , 

Q</3 

where laf3 := \za — zpl- It is easy to see that equations © yield (0) and 
hence the system has a Hamiltonian structure. In particular H (the total 
energy) is a conserved quantity. From it follows that also 

Z = ^aZa 
a 

and 

I = '^'^pI^ Pa = \Za\ 
a 

are conserved quantities. 
^Cf. §1.1.4] 
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Let r = Ta- If r / we call Z/T the barycenter of the system. A 
calculation shows that 



M:= ^r,r^/2^ = r/-|z 



a<l3 



Finally, from it also follows that the angular momentum or virial 



is an invariant of the system, i.e., a constant depending only on the vortici- 
ties. Without loss of generality we will assume that Z = when T ^ 0, and 
hence M = TI. 

Notice that all initial configurations that end in collision must satisfy 
M = I = (even when F = 0). This can be viewed as an extension to the 
three vortex problem of Sundman's theorem from celestial mechanics, which 
states that total collision implies that the angular momentum must be zero, 
see (Tj. 

We say that the dynamical evolution of N vortices on the plane is a self- 
similar trajectory if it can be expressed as Zi{t) = ("(t)zj(O), i = 1, . . . N, 
for some C : ^ — > C (hence Z = 0). As we will see below, a self-similar 
trajectory that is not a relative equilibrium (defined at the end of section I^J 
leads to total collision or ejection. 

2 Reduced coordinates and trajectories leading to 
total collision 

A given configuration of N point- vortices can be recovered if one knows 
all the distances lap between the vortices and the overall orientation of the 
system. The reduced system obtained after ignoring the orientation is called 
the shape configuration. This will be enough for studying candidate 
collision trajectories. 

From now on we will concentrate on the case = 3. The way we 
prove the main result of this section, stating that any collision of planar 3- 
point- vortices is always self-similar, is by writing the conditions of collapse 
in terms of coordinates introduced by Novikov |Hj, which are essentially 
mutual distances, without requiring that the virial (or angular momentum) 
vanishes. We prove that the cases of positive or negative virial leads to no 
solutions where all the 3 vortices approach the origin. For this we apply a 




a 



a<l3 
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projective blow-up on a plane in the space of parameters. The only case left 
is that of null virial, which is known to contain the solutions of self-similar 
collapse. 

Let C3 denote the cyclic permutations of (1,2,3) and let us introduce 
the notation 6j := {i,j,k) £ C3. 

The shape configuration of the system is a point (61,625^3) in the first 
octant of such that the triangle inequalities 

V^+7^>\/6fc, {i,j,k)eC3 (3) 
are satisfied. These three inequalities are equivalent to the single inequality, 

6? + 6^ + 6i < 2(6162 + 6163 + 6263) • (4) 

Geometrically, the set of points (61,62,63) G satisfying 6j > and © is 
a solid cone which we call the admissible cone, with symmetry axis the 
line 61 = 62 = 63; this is clear from the symmetry and homogeneity of @. 

Now we study the trajectories determined by the constraint that H and 
M must be conserved quantities. As it was mentioned in the introduction, 
all candidate trajectories for total collision must satisfy M = 0. Another 
condition for total collision is that not all the Fj have the same sign; for if 
that were the case then H would tend to 00 as lai3 0. For concreteness 
let us assume that F1F2 > 0, F1F3 < and thus F2F3 < 0. 

Let h{t) = (6i(t), 62(t), 63(t)) be a dynamical trajectory of the system 
satisfying Since H and M are conserved, we have that, for all t, 

^ = --3-(FiF2ln63 + F2F3ln6i + F3Filn62) , 
47r 

= F1F263 + F2F361 + F3F162 . 
Exponentiating the first equation we obtain the equivalent system 

63 = /i6™6^, (5a) 
m6i + n62 — 63 = (5b) 

where m = -F3/F1, n = -F3/F2, h = e'^^'^/^^^'^^'i > 0. Equation ^ 
expresses conservation of energy and equation H5b|) the condition M = 0. 
Notice that m, n and h are strictly positive. 

The dynamics of the system is contained in the region T = CD {M = 0}, 
where C is the admissible cone determined by ^ and {M = 0} denotes the 
plane determined by H5b|) . If T 7^ then T is a planar region contained in 
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the first octant of bounded by two rays ri, r2 through the origin (possibly 
ri = r2). The projection of these rays on the x-y plane determine two rays 
with slopes pi, p2 such that < pi < P2 1^ oo. Its values are obtained by 
eliminating 63 from the equality case in Q) and equation ()5bl) . which gives 
the quadratic equation in p = 62/^1 ) 

(n — + 2{mn — m — n — l)p + (m — 1)^ = , (6) 

whose roots are pi^2 = (1 =t y^)^/(n — 1)^, where 

/? := m + n — mn . 

For /3 < 0, /3 = and /3 > 0, with /? 7^ 1, equation @ has no real roots, 
exactly one positive root and exactly two positive roots, respectively. Notice 
that /3 = liffm = lorn = l, which corresponds to © having a root p = 
or p = 00, respectively. Thus, T is not empty iff /3 > 0, it degenerates to a 
line when /3 = and has a boundary line on the cartesian plane x-z or y-z 
iff m = 1 or n = 1, respectively. 
Let 

, = l-n.-n = r3(i^ + ^ + ^) . (7) 
We now consider the two possible cases: 



Case 7 = 0. In this case V = 0, which corresponds to the well known 
self-similar collisions described in Notice that in this case < (3 < 1, 
and that T contains the symmetry axis of C. From the homogeneity of 
it is clear that for each ray r through the origin contained in T there is an 
energy value h such that r is a set of solutions of ©• Therefore, the region 
T is foliated by rays through the origin which are the dynamical orbits of 
the system, except for the boundary rays ri, r2 and the bisector of T, which 
consist of equilibrium points, as we show in section |1J 

Case 7/0. Letting x = 61, y = 62, z = 63 we obtain from equations © 
the equivalent conditions 

(m + np)xi-™-" - /ip" = , 

(8) 

y = px , z = mx + ny . 

Therefore, 

h}-hp^/i /ii/7p{i-m)/7 h^hp^h 

X v = z (9) 

(m + np)i/7 ' ^ (m + np)V7 ' (m + np)(™+")/T ' 
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Equations © describe a parametrized curve c{p) = {x{p),y{p), z{p) defined 
and different from zero for all p £ (0, oo). If 7 > then (x, y,z) ^ as p — > 
or p ^ 00 (hence the curve is bounded). If 7 < we see that, if p — > 0, 
respectively p ^ oo, then z behaves asymptotically as p"'''"', respectively 
p-m/'r^ In either case z^ooasp^Oorp^oo (hence the curve is 
unbounded). 

The dynamical orbit is contained in the portion of c lying on the planar 
region T. As we will show in section |11 the only equilibrium points occur 
when /3 = and in this case T degenerates to a line of equilibrium points. 
If /3 / and /3 / 1 then, from the discussion following equation the 
dynamical orbit consists of the curve given by Q with < pi < p < P2 < 00. 
This is a simple curve away from the origin joining two points, one on each 
boundary ray of T, with the property that each ray through the origin in T 
intersects the curve exactly once. If /? = 1 then, since in this case m = 1 or 
n = 1, 7 < and the curve does not tend to the origin as p — > or p ^ cxd. 

We have thus proved: 

Theorem 1. For a system of three point-vortices on the plane all total 
collisions are self-similar. They occur only when vorticities satisfy V = 
and they lie on the level M = 0. 

Remark. The parametrization of the dynamical orbit by the variable p 
can be interpreted as a projective blow-up in the x-y plane. 

Remark. The equilibrium configurations 61 = 62 = ^3 = correspond to 
rigid motions of the vortices (all mutual distances are preserved). It is usual 
to call these motions relative equilibria. 

3 Regularization of collinear configurations 

In this section we introduce new coordinates that regularize the vector field 
at the boundary of the cone C given by @. Also, we show that it is not 
possible that two vortices collide while the third one remains away from 
the collision. (Hence, the only kind of collision for three planar vortices is 
self-similar collapse.) 

First we obtain^ the evolution equations for the 6j's. Let wt = zj — z^, 
hence bi = Here and in what follows (a, /3,7) is a cyclic permutation 

^See |5] for a geometric derivation of these evolution equations. 
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of (1,2,3). Prom the equations of motion we get 



-. = f (^ + ^-^ + ^) (10) 
zvr \ Wk Wi Wj Wi 



Thus, 



d , , . 

— [WiWi ) = WiWi + WiW. 



Im — + 

vr \wj WkJ 



Noticing that Wi = —{wj + Wk), 



vr \ WjWk J \Wk Wj 

1 1 



= — Im {WjWk, 1 -. - -. 

vr \wk\ 
In the last equality we have used the identity 

Im ( — + — ) = lm{ziZ2) ( -tAo - -rKn 
\Z2 zij \\ziV F2r 

for all zi, Z2 G C \ {0}, which is easily proved by expressing the Zi in polar 
coordinates. Let A denote the oriented area of the triangle formed by the 
point vortices located at zi, -^3 £ C. It is easy to see that 2 A = lm.WjWk- 
Hence, 

6. = ^(i-l) , J, A:) cyclic. (11) 

This represents in fact two vector fields in the interior of C, since A is 
defined up to a sign. It becomes zero at the boundary of C, but in a singular 
way since the expression of A in terms of the 6j's involves the square root of 
a function equal to zero at dC (cf. equation (|n|) below). Thus, even though 
the original equations of motion ^ are regular around generic collinear 
configurations, the equations of motion (|llj) on the reduced space of shape 
configurations are not. 

The vector field defined by becomes infinite when 6j — > while hj , hk 
remain away from zero, for some (i, j, k) G C3. But we now show that such a 
binary collision can not be reached. Indeed, assume that 6j — > as t — > toj 
to < 00. From the triangle inequalities Q we have that y/hl > \ -\/hj — y/bk\ 
and hence \bj — hk\ — > 0. Expressing the conservation of energy as 

h = b]'WJWf^ (12) 
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we see that bj,bk can not tend to c with < c < oo. Therefore either 
bj,bk ^ or bj,bk oo, as t — > io- The former corresponds to total 
cohision, the kind considered in the previous section, and by definition is 
not a binary colhsion; so let us assume the latter. Conservation of M can 
be written as 

M = TiiTj + Tk)bk + TjVkh + TkT.ibj - bu) . 

Since the last two terms in the right-hand-side tend to zero while b^ ^ oo 
we conclude that M = and Tj + Tk = 0- Therefore ((T^ can be rewritten 
as 



bk 

The left-hand-side tends to or oo, while the right-hand-side tends to /i > 0, 
which is a contradiction. We have thus proved: 

Theorem 2. The dynamical trajectory of any three-point-vortex configura- 
tion does not evolve into a binary collision. 

Now we introduce coordinates (a, A, 6) which, as we will show, regularize 
the vector field at the boundary of the admissible cone The price that 
we will pay is that in the new coordinates the vector field at the equilateral 
configurations (which, as we show in section |1J are equilibria) will become 
undefined. These coordinates are defined as follows. Let 

A = 61 + 62 + ^3, 6j = |zj - ZfcP ■ (13) 

Consider three points on the plane (|T3|) for a fixed A: P2 = (A/2, A/2,0), 
Pi = ('^/S, A/3, A/3) and b = (61, 62, 63). Thus Pi is the center of the disc V 

on the plane (|T3|) bounded by the cone (jU and P1P2 is the radial segment 
of V perpendicular to the 6162-plane. Let 6 be the angle needed to take the 

ray P1P2 to Pib. Finally, let a = 12A/(\/3A). 

From Heron's formula for the area of a triangle. 



4A = eJ2{bib2 + bibs + 6263) -bj-bl-bj 

(14) 



= e^A2- 2(62 + 61 + 61), 

where e = ±1, depending on the orientation of the triangle formed by the 
point-vortices. It is clear that \A\ takes its maximum value when 61 = 62 = 
63 = A/3 so that l^lmax = \/3A/12. Let 

M := (-1, 1) xR+ X S'^ = M-UMoU M+ 
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where 

M- := (-1, 0) X M+ X 5^ M+ := (0, 1) x M+ x 

and 

Mo := {0} X E+ X 5^ 

is the common boundary of 7W_ and A^+. It is easy to see that the cor- 
respondence between the variables (a, A,^) and {bi,h2-,b^) induces a difFeo- 
morphism between M± and C \ £, where C is the admissible cone defined 
by and 

£ := {(61,^2 5^3) G C I bi = 62 = ^3} (equilateral configurations) , 

and a diffeomorphism between Mq and dC\£ (these are the collinear config- 
urations minus total collapse). In contrast, the boundary dAi corresponds, 
in a non-bijective manner, to the equilateral configurations £. (Notice that 
the total collision belongs to £.) 

Let us obtain the expressions of the bi as functions of a, A, 9. From (|14l) 
and the definition of a we have that 

bl + bl + bl = ^{3-a'). (15) 

Also, 

Pib ■ P^2 bi + b2- 263 



cos I 



pTb 1 1 pTp2 I VeVb'i + b'i + bi-xys 

Using that sin^ = 1 — cos^O and that sign(sin 6) = sign(62 — 61) we get that 

62 - bi 



sm ( 



^/2./b^+h + h^W/?, ' 
Using (|15() these expressions can be simplified: 

A -363 . . \/3(b2 - 61) 

cost^ = -; — , smb* = — ; — ^ . (16) 



A^/r^T? A^/T 



From (|13|) and (|16j) we obtain a linear system for the 6j's whose solution is 
A 



bi = - ( 2 + a/i - a2(cos6' - \/3sin( 
6 V 



'l- 




'l- 





62 = -^(2+ Vl-a2(cos0 + V3sin0)) , (17) 
A 



63 = 3 (l- \/l-a2 cos^j , 



10 



so that the ftj's are differentiable functions of (a, A,^) G M.. 

From (fTT]) . l(T^ and the definition of a we see that A = a fx, where fx is 
a rational function of the ftj's, with parameters Fj. From and (|14() we 
get that 



A 



1 



Fi ( 1 - 1) (-61 + 6, + 63) + ^2 - ^ ) (61 - 62 + hj 



+ r3 



1 1 

61 62 



(fei + 62 



• (18) 



Hence, a = (12/^/\/3 — a'^fx)/X, where is another rational function of 
bi,b2, 63 with parameters Fj. Both /a and f^ are differentiable except when 
bi = for some i. Hence the same is true for A, q and, because of and 
H17() . it is also true for 9 = cos9{sm9)' — sm9{cos9)' . 

Therefore the variables (a. A, ^) regularize the vector field (lllj) on M := 
M\B, where 

e = {(0, A, 2n7r/3) G 7W | n = 0, 1, 2} 

is the set of binary collisions, i.e. bi = for some i, which is a codimension 
2 closed subset of A^. In this set the vector field is singular and, as stated 
in theorem 12 it can not be reached even asymptotically. 

The vector field on Ai can be easily computed, as we now show. From 
(jlTp one computes the total derivative D{bi,b2,b3)/D{a, X,9), whose deter- 
minant is equal to y/SaX'^ /18, so that it is invertible on 7W iff a 7^ 0. Its 
inverse is computed to be 



B{a,X,9) VT 



0(61,62,^3) 



aX 



B 



where 



/VT- 



B 



Vl-a^-v^a 



cos I 



\ 



-Q!(\/3 sin 9+3 cos i 



-a(\/3 sin 6—3 cos I 



V3(l-a2) 

Thus, from and the chain-rule. 



V3(l-a2) 



2a sin 9 
l-a2 



a" 



247r 







B F2(63 1 




\r3(6r^ 
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where the bi are given by lfT7|) . As expected, the vector field {a,X,6) is 
regular on A^. In particular, we have achieved the regularization of the 
vector field around the collinear configurations (excluding the non dynamical 
binary collisions). 

Remark. It is clear that A = iS a = 0. Noticing that from = 
^^^2^^ we see that when ^ = 0, we also have that A = iff d = 0. 

4 Equilibrium configurations 

In this section we determine the equilibrium configurations of the system de- 
fined by equations This completes the discussion in section |21 where the 
description of the dynamical orbits requires distinguishing the equilibrium 
points. As remarked in that section, these equilibria have the interpretation 
of being the relative equilibria of the 3-point vortex system. 

Directly from we see that a necessary condition for equilibrium is 
that either 6i = 62 = &3 (equilateral configurations) 01 A = (collinear 
configurations). Since the vector field (|11() is smooth around equilateral 
configurations, apart from total collision, then being at an equilateral con- 
figuration is also a sufficient condition for equilibrium. On the other hand, 
A = is not a sufficient condition for equilibrium because the expression 
of A in terms of the ftj's involves a square root and hence the vector field 
is not smooth around collinear configurations. In this case the equilibrium 
conditions at those configurations can be easily determined with the help 
of the regularized vector field (d. A, 6) of section |3| from the expressions for 
(d, X,9), B and the remark at the end of section|31we see that if yl = then 
a) A = ^ = and b) d = iff i = 0. Therefore, 

Lemma 3. The equilibrium configurations of the evolution equations (|1H) 
are a) the equilateral configurations and b) the collinear configurations for 
which ^4 = 0. 

We now determine the conditions for ^ = when A = in the case 
when M = 0, which is the context of section [2 First we observe that, from 
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H5b|l and ((TH|) . A can be expressed in terms of p = 62/61 by 

cA = —6123 + -6231 - ^312 , 
m n 

h , h -1 , ^2 

— = m + np , — = mp + n , -z— = p , 

Ol 62 61 

where c = — 327r/r3 and k) G C3. The values of p that correspond to 
^ = (the boundary hues of T) are given by the roots of © • A calculation 
shows that substituting these values of p in ()19() and expressing m and n in 
terms of /3 = m + n — mn and 7 = 1 — m — n yields 

Therefore i = iff either /? = or 7 = 0. Noticing that /3 = -rTr^r and 
recalling, from ((7|, that 7 is proportional to the virial, we conclude: 

Proposition 4. Assuming M = 0, a collinear configuration is an equilib- 
rium if and only if either the total vorticity or the virial vanishes. 

Remark. From the discussion following equation @ we see that if the to- 
tal vorticity vanishes then T, the set of states satisfying M = 0, degenerates 
to a single line of equilibrium points. 



We conclude by noticing that the self-similar trajectories of the case 
M = 7 = discussed in section |21 correspond to collisions or ejections 
depending on the sign of A, where A is one of the regularizing coordinates 
introduced in section |3| From the definition of A and equation we 
compute that 

27r63 ; _ ((2 - m)p + m + 1) (p - 1) 



where p = 62/61. The right hand side is zero iff p = (m + l)/(m — 2) or 
p = 1. The former can not happen since it is incompatible with p > and 
< 771 < 1. Also, for these ranges of p and m, (2 — m)p + m + 1 > 0. There- 
fore, if ^ / 0, A = iff p = 1; and if p 7^ 1 then sign(A) = sign[r3 A (p—l)] = 
sign[r3 a sinO]. Thus, for a given orientation of the vorticities, the bisector 
of T (which corresponds to the equilateral configurations) separates the col- 
lapsing and expanding self-similar trajectories. In terms of the regularizing 
coordinates, we see that the collapsing or expanding nature of the self-similar 
trajectories changes when a —a or 1— > —9. 
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